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Abstract 

This  report  is  a  summary  of  a  two-year  effort  by  the  PFs  research  team  and  collaborators.  We  have 
demonstrated  the  prediction  capability  of  module-based  multiscale  ionized  gas  (MIG)  flow  finite-element 
code  for  mitigating  inert  gas  flow  separation  using  rf-driven  dielectric  barrier  discharge.  We  extended  our 
work  towards  developing  a  multidimensional  first  principles  theoretical  model  of  the  non-equilibrium  real 
gas  discharge.  The  electric  force  field  generated  by  asymmetrically  arranged  plasma  actuators  operating  at 
rf  is  self-eonsistenlly  coupled  to  the  gas  dynamics  of  the  air-like  N2/O2  gas  mixtures.  Specifically,  several 
physical  and  geometric  parameters  such  as  the  amplitude  and  shape  of  excitation,  dielectric  constants, 
initial  ionization  level,  and  electrode  shape  are  studied.  The  effect  of  polyphase  power  supply  to  an  array 
of  actuators  has  also  been  explored.  It  is  found  that  (a)  favorable  ranges  of  frequency  arid  electric  Held 
exist  for  most  effective  flow  control,  (b)  the  momentum  transfer  to  the  neighboring  gas  is  cumulative  for 
an  optimum  phase  angle  and  distance  between  the  exposed  electrodes.  Also,  stream  wise  momentum 
transfer  does  not  keep  on  increasing  with  the  number  of  actuators.  We  have  worked  in  close  collaboration 
with  Dr.  Datta  Gaitonde  (VAAC/AFRL)  to  integrate  MIG  with  FDL3D1  code  and  to  simulate 
electrodynamic  mitigation  of  three-dimensional  wing  stall  about  a  symmetric  airfoil.  Similar  collaboration 
with  VAAC/AFRL  researchers  is  underway  towards  modeling  plasma  control  of  transitional  vortical  flow 
over  a  low  sweep  delta  wing.  These  efforts  set  the  basis  for  applying  the  simulation  method  for  actively 
controlling  vortical  flows  past  common  aerodynamic  configurations  like  a  vehicle  forebody. 
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I.  Introduction 


Recent  successful  experiments  and  simplified  numerical  models  for  boundary  layer  flow  actuation  using 
dielectric  surface  discharge  at  atmospheric  pressures  show  tremendous  potential  for  aerospace  and  many 
other  applications  [1-6],  These  watt-level  actuators  can  efficiently  apply  large  electrohydrodynamic 
(EHD)  forces  in  a  relatively  precise  and  self-limiting  manner,  have  rapid  switch-on/off  capabilities  and 
require  no  moving  parts.  A  detailed  theoretical  understanding  of  the  actuation  mechanism  is  necessary  to 
design  these  devices  for  more  effective  momentum  injection  for  better  flow  control  Typical  electrode 
arrangements  for  these  actuators  have  a  grounded  electrode  which  is  surrounded  by  an  insulator  and  a 
voltage  fluctuating  with  radio  frequency  (rf)  is  applied  to  the  electrode  exposed  to  the  gas*  The 
experimental  data  suggest  1-20  kV  peak-to-peak  applied  voltage  with  2-50  kHz  rf  is  suitable  for  these 
actuators  operating  at  atmospheric  pressure.  The  plasma  discharge  at  this  pressure  is  highly  collisional 
causing  an  efficient  energy  exchange  between  charged  and  neutral  species.  Numerical  simulation 
predicted  new  insight  into  the  near  wall  momentum  exchange  physics.  It  is  observed  that  dielectric 
surface  becomes  negatively  charged  for  most  of  the  power  cycle  and  a  time  averaged  force  acts  on  the 
plasma  predominantly  downstream,  with  a  transverse  component  towards  the  wall  Momentum  of  the 
plasma  couples  to  real  gas  through  collisions,  which  enhances  near-wall  momentum  yielding  a  wall-jet¬ 
like  feature  for  quiescent  gas.  A  suitable  range  for  electrical  and  physical  control  parameters  has  also  been 
identified  for  which,  the  induced  momentum  most  effectively  mitigates  separated  flow  over  a  flat  plate 
and  the  wing  stall  characteristics*  It  is  found  that  the  momentum  transfer  to  the  neighboring  gas  is 
cumulative  for  an  array  of  actuators  showing,  however,  a  diminishing  return  with  increasing  number  of 
electrodes. 

M  Objective 

This  numerical  simulation  effort  aims  to  explore  the  transient  piasma-gas  interaction  processes  of  a  near 
electrode  RF  excited  dielectric  barrier  (optical  glow)  discharge  for  helium  gas  and  air-like  N i/O?  mixture 
at  atmospheric  conditions.  These  processes  have  been  shown  experimentally  to  have  profound  impact  on 
both  high-  and  low-speed  flows  relevant  to  the  Air  Force  interest.  The  research  goal  is  to  develop  from 
first  principles  a  self-consistent  theoretical  model  of  suitable  complexity,  to  implement  it  into  a  powerful 
robust  and  efficient  numerical  tool,  and  finally  to  demonstrate  the  effect  of  surface  plasma  discharge  on 
real  gas  flow  aboui  a  realistic  geometry.  Comprehensive  understanding  of  the  pertinent  mechanisms  will 
significantly  impact  the  high  speed  flight  endeavor;  enable  introduction  of  new  concepts  exploiting 
fundamental  EHD  flow  control  and  guide  current  and  future  experimental  efforts. 
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111.  Approach 


Due  to  a  combination  of  electrodynamic  and  collisional  processes  the  proper  mechanism  of  which  is  yet 
to  be  theoretically  resolved,  charge  separated  particles  induce  the  gas  particles  to  move.  For  accuracy  and 
fidelity,  it  is  essential  that  the  force  relations  be  derived  from  first  principles  through  a  simulation  of  the 
elementary  mechanisms  that  yield  the  discharge.  We  are  developing  [4-9]  such  models  in  the  following 
manner.  First,  a  drift  diffusion  partially  ionized  gas  model  for  atmospheric  air-like  N2/02  gas  mixture  is 
formulated  for  nine  species.  Various  reaction  rate  coefficients  [10]  are  obtained  based  on  electron 
temperature  calculated  from  the  instantaneous  electric  field.  The  reactions  include  momentum  transfer, 
molecular  ionization,  dissociation,  dissociative  recombination  and  secondary  electron  emission. 
Atmospheric  ratio  of  3.6  is  taken  for  nitrogen  to  oxygen  gas  molecules.  The  model  is  anchored  in  the 
module  based  multiscale  ionized  gas  (MIG)  flow  code  [4-6,  11-13].  The  self-consistent  formulation  is 
solved  using  a  Galerkin  variational  formulation  based  finite-element  method  [12,13]  to  obtain  electron, 
positive  and  negative  ions  and  neutral  species  densities  of  nitrogen  and  oxygen,  and  the  electric  potential 
in  two  spatial  dimensions. 

The  BHD  body  force  thus  calculated  from  separated  space  charge  and  electric  Held  is  then  employed  for 
solving  gas  hydrodynamic  equations.  Barrier  discharges  in  air-like  N2/02  mixtures  consist  of  a  multitude 
of  micro  and  nanoscale  discharges  [8f  MIG  implements  a  powerful  high-fidelity  finite-element  (FE) 
procedure  adapted  from  fluid  dynamics  to  overcome  the  stiffness  of  the  equations  generated  by  multi- 
spec  ics  charge  separation  phenomena.  Implemented  FE  techniques  are  especially  suitable  for  operating  at 
disparate  spatial  and  temporal  scales  and  for  their  adaptability  to  arbitrary  multidimensional  geometries 
and  boundary  conditions  [13].  Here,  a  2D  biquadratic  basis  polynomial  is  employed  to  interpolate  the 
state  variables  along  with  implicit  time  marching  to  solve  the  continuity  and  momentum  equations  for  all 
species  and  gas.  The  solution  process  consists  of  two  steps  -  first  solve  for  all  species  (including  neutrals) 
lo  calculate  the  charge  separation  and  electric  field  up  to  a  periodic  asymptote  with  a  timestep  of  up  to 
100  times  the  dielectric  relaxation  timescale,  then  use  those  as  body  force  lo  predict  the  air  velocity  and 
continuity  evolving  in  time  with  the  flow  timescale.  The  solution  of  the  Newton  iteration  for  the  resulting 
matrices  is  convergent  at  any  given  time  step  when  the  maximum  value  of  the  residual  relative  L?  norm 
for  each  of  the  state  variable  becomes  smaller  than  a  convergence  criterion. 

Several  cases  including  various  electrode  configurations  and  polyphase  power  input  are  simulated  using 
both  helium  and  air  as  working  gas.  The  solutions  computed  are  then  qualitatively  benchmarked  with  the 
reported  test  data.  The  computed  END  force  data  are  also  transitioned  to  the  VAAC/AFRL  researchers 
who  predicted  the  flow  control  mechanism  over  a  stalled  NACA  0015  wing. 
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IV.  Discharge  Formulation 

The  following  two-dimensional  three  species  collisional  plasma-sheath  model  is  solved. 


Charge  continuity: 

dr  dx, 

(la) 

Charge  momentum: 

(lb) 

Potential: 

(2) 

Neutral  continuity: 

On  dnKVni 

— ”  +  — 2_2.  =  _n  -  +  rn  n 

dr  dx, 

(3) 

The  charge  particle  a  =  ej  distributions  are  considered  non-Maxwellian  and  inertia  terms  are 
neglected.  The  electron  temperature  is  nearly  uniform  at  leV  =  ll,600K  and  the  ions  and 
neutrals  are  in  thermal  equilibrium  at  300K.  The  working  gas,  helium,  is  maintained  at  bulk 
pressure  p  =  300  lorr  and  temperature  T=  300K. 

Transport  properties  are  taken  from  the  literature.  The  electron  diffusion  is  obtained 

from  Einstein  relation,  D-  where  Te  is  the  energy  in  electron  volts,  e  is  the  elementary 

charge,  zris  the  permittivity,  and  pe=e/(m^veh)  is  mobility  of  an  electron,  where  i^adO'Ts  is  the 
electron-neutral  collision  frequency.  We  emphasize  that  transport  properties  used  from  the 
published  literature  are  to  be  taken  as  nominal  values  facilitating  the  development  of  the 
numerical  framework,  rather  than  providing  new  thermo-chemical  data  for  these  discharges. 
Similarly,  the  ion  diffusion  D,=500  cm2/s  at  300K.  and  the  ion  mobility  //,  is  given  as: 


pp,  =8xl01(l-8xl0‘1£/ p)  cnTV's  'lorr,  for  El  p<  25  V  cm  'ton- 1 


PH,  = 


4.|x!0J 


JFTp 


I- 


27.44 


{E*P) 


cm2 V'Vtorr,  for  E!  p>  25  V  cm  'torr  ' 


(4) 
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where  8  ^ '  +  is  the  electric  field  magnitude  with  components  ~  an(j  p  js  the  gas 

pressure.  The  index  j  takes  the  values  x  and  y.  The  ionization  rate  z  for  helium  gas  used  here  is: 

/ 

—  ppcE  s4;  A  =  4.4  cm  ‘torr4  and  B  =  14  [V/(cm  ton*)]04 

(5) 


z  =  A  exp 


[WpT 


where,  ne  is  the  electron  number  density.  The  coefficient  of  recombination  on  the  right  hand  side 
of  Rqn.  (la)  and  (3)  is  given  as,  r  =  {yM,<rr„{VlA))  =  1.09x|0~1°7'",'\mVs.  Here,  Ye\b  is  the  electron 

thermal  velocity  and  crei  is  the  electron-ion  collision  cross-section,  for  simplicity,  we  have 
presently  ignored  recombination  effect  on  momentum,  as  well  as  secondary  emission,  at  this 
stage. 


The  dynamics  of  neutral  particles  are  determined  using  the  following  neutral  momentum 
equations: 


8V*>.  +  v  dV”'  -  m  d 


dt  dx,  mnnn  dxt 


(  AA  \ 


dx,  J 


+  ^sHk. v  (j V  -V  )  +  -^-v  (V  -  V  )  +  n  zV  In  (6) 

trrt  v  i?J  "/  /  w  \  tj  m  }  €  ill  n  v  * 


The  averaged  velocity  of  the  iluid  of  density 
momentum  equation: 


itv  P~mSn„+n.)  \ 


is  based  on  the  following 


dVe  „  dV  m  r  T  op  8 

“  +  —  =  —  e{»,  -K)E,  — — -T-  +  —  i 
oi  cx}  p  p  dx  f  dxf 


(7) 


where  t  ,  =  —  C-^-  +  ^  —  with  p  as  the  gas  viscosity  (=2.066x  1 0°  Ns/m" )  and  6,/  is 


\ 


\  J 


3  p  dxk 

the  Kronecker  delta.  The  factor  tn  introduced  to  modulate  the  effect  of  etectric  Held  is  set  to 
unity. 


The  multiscale  ionized  gas  (MIG)  flow  code  was  used  to  solve  Eqns.  (l)-(7).  MIG 
implements  a  versatile  high-fidelity  finite-element  (FE)  procedure  adapted  from  fluid  dynamics 
to  overcome  the  stiffness  of  the  equations  generated  by  multi-species  charge  separation 
phenomena,  FE  techniques  are  especially  suitable  for  their  adaptability  to  arbitrary 
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multidimensional  geometries  and  boundary  conditions.  Mere,  a  2D  bilinear  finite  dement 
formulation  is  employed  together  with  4  order  Runge-Kutta  time  marching  to  solve  the 
continuity  and  momentum  equations  for  all  species  and  fluid.  The  solution  process  consists  of 
two  steps.  First,  using  Eqns.  (la),  (2),  (3)  and  (6),  a  global  matrix  is  formed  and  solved 
simultaneously  obviating  the  need  for  any  special  sub-iteration  for  the  Poisson  solver.  The 
species  density  and  charge  velocity  thus  calculated  are  then  used  in  Eqn.  (7)  to  predict  the  fluid 
velocity.  The  Galerkin  weak  statement  associated  with  a  variational  integral  underlines  the 
development  of  this  numerical  algorithm.  The  solution  of  the  Newton  iteration  is  converged  at 
any  given  timestep  when  the  maximum  value  of  the  residual,  relative  norm  for  each  of  the  state 
variables,  becomes  smaller  than  a  chosen  convergence  criterion  of  1  O'". 

The  computational  domain  (x,y:  0,0.01m)  consists  of  a  lower  half  (y :0,0.005m)  insulator  with 
dielectric  constant  e<j=  4.5  eq,  and  a  upper  half  filled  with  inert  helium  gas  of  Ef  =  1.0055  Eo 
where  Eo  is  permittivity  of  vacuum.  Inside  the  insulator  the  current  due  to  motion  of  charged 
particles  is  forced  to  zero  while  the  displacement  current  is  balanced  with  the  total  current  at  the 
gas-dielectric  boundary.  The  schematic  in  Figure  I  shows  two  electrodes  in  which  the  bottom 
electrode  is  grounded.  $0)=  0,  and  a  sinusoidal  RF  alternating  potential  <prf=tf>rms$\n27ift  with 
^,ns=l  kV  and f=5  kilz  is  applied  to  the  lop  electrode.  Each  electrode  is  infinitesimally  thin  and 
2mm  long. 

In  order  to  identify  the  effect  of  geometric  configuration  two  different  electrode  arrangements 
were  simulated.  In  the  first,  the  two  electrodes  were  kept  vertically  aligned.  Hereafter,  we  will 
denote  this  as  the  symmetric  configuration.  In  the  second  case,  the  electrodes  arc  offset 
horizontally  by  2mm.  I  bis  will  be  noted  as  the  asymmetric  configuration.  A  decrease  in  the 
offset  horizontally  and/or  vertically  increases  the  electric  field  for  the  same  applied  rms  potential 
which  in  effect  increases  induced  electric  force.  For  all  simulations  electrons  are  assumed  to  be 
isothermal  at  the  boundary  and  maintained  at  leV  (-1 1.600K)  while  the  ions  are  cold  (300K)  at 
300  lorr.  The  electron  flux  at  the  exposed  electrode  is  based  on  the  electron  thermal  velocity  and 
is  directed  towards  the  wall.  Homogeneous  Neumann  boundary  condition  ( dn/dx=0 )  is  imposed 
for  ions  at  the  exposed  electrode.  The  solutions  are  verified  by  qualitative  comparison  with 
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published  results.  The  simulations  are  then  employed  to  explore  the  enhancement  of  near  wall 
fluid  velocity. 

The  ionized  helium  gas  is  numerically  modeled  using  the  finite-element  based  MIG  flow  code. 
The  code  is  modular  and  separate  subroutines  can  be  written  to  model  different  physics.  Here, 
the  equation  sets  (1)  (6)  can  be  written  with  operator  L  as  L  (q)  =  0  where  q  contains  state 

variables  like  number  densities,  velocities  and  potential.  Multiplying  with  a  permissible  lest 
function  //and  integrating  over  the  spatially  discretized  domain  Q,  the  variational  statement 
results  in  the  weak  form 


WSh  =  3. 


jl  rjL(q)dr] 


=  0 


for  a  discretization  h  of  f7  =  Uf7,  and  is  the  non-overlapping  sum  over  the  elements.  Thus,  for 
example,  the  GWS  form  of  Eq.  (2)  becomes. 


•  \t\  n,  a 


drf 


(8) 


where  F+  is  the  solution  residual. 

The  Jacobian  matrix  J-  [dF/8q\  in  the  global  [,/j.{a?}  =  -{F)  is  resolved  using  LU-decomposition 

scheme  for  updating  change  in  discretized  solution  vector  q  at  each  iteration.  The  terminal  non¬ 
linear  ordinary  differential  equation  (ODE)  systems  are  solved  using  implicit  Euler  method  for 
temporal  evolution  and  N-R  iterative  algorithm  for  the  non-linear  matrix  algebra.  The 
convergence  criterion  for  all  variables  at  any  iteration  is  set  at  ID'4.  Solution  stability  is  ensured 
by  appropriate  selection  of  time  marching  step  size  and  the  introduction  of  artificial  diffusion. 
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V.  Major  Developments 


I.  Helium  discharge  model  for  quiescent  How 

For  6  kHz  and  1.6  kV  peak-to-peak  with  alumina  insulator,  model  predictions  for  charge 
densities,  the  electric  field  and  gas  velocity  distributions  are  shown  to  mimic  trends  reported  in 
the  experimental  literature.  We  also  predict  the  electron  charge  accumulation  on  the  dielectric 
surface  (Figure  la.b)  self- limiting  the  discharge  and  influencing  the  essentially  unidirectional 
charge-neutral  momentum  exchange  forming  a  wall  jet  (Figure  1c).  The  velocity  field  based  on  a 
quiescent  flow  initial  condition  shows  stream  wise  component  of  the  computed  gas  velocity  in 
Figure  Ic  and  depicts  a  strong  wall  jet  downstream  of  the  rf  electrode  away  from  the  dielectric 
surface.  The  highest  momentum  transfer  between  the  gas  and  the  charge  particles  occur  slightly 
downstream  of  the  peak  force  location  shown  in  Figure  Id.  The  local  vertical  line  plots  at  six 
streamwise  locations  describe  how  the  flow  velocity  increases  responding  to  the  increase  in  axial 
electric  field  and  experimental  results  show  very  similar  feature  [I].  The  peak  jet  velocity  is  -1.7 
times  the  calculated  freestream  velocity. 


(a) 


Figure  I :  Plasma  actuation  of  quiescent  helium  gas.  (a)  Charge  particle  distribution  shows 
electron  accumulation  on  the  dielectric,  (b)  time  average  of  charge  density,  (c)  streamwise  gas 
velocity  profiles  at  different  locations  along  the  flow,  (d)  force  variation  at  the  peak  of  the  cycle. 


2.  Helium  discharge  model  for  separated  flow 

Higher-fidelity  models  are  included  to  yield  a  more  sophisticated  framework  to  predict 
discharge- induced  momentum  exchange.  Here,  the  complete  problem  of  a  dielectric  barrier 
discharge  based  separation  control  with  axially  displaced  electrodes  is  simulated  in  a  self- 
consistent  manner.  The  body  force  is  computed  self-consistently  from  the  three-body  electro  gas 
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dynamics  using  the  applied  rf  potential  and  the  dielectric  properties  of  insulator  and  gas  media 
based  on  first  principles.  For  5  kHz  and  l  kV  with  Kaplon  dielectric,  specific  details  of  the 
computed  flow  components  and  streamwise  force  is  given  in  Figures  2a-2c.  While  the  normal 
force  component  is  primarily  downward,  the  time  average  of  the  streamwise  force  component  in 
Figure  2d  shows  the  dominance  of  the  positive  force  to  push  the  fluid  in  the  forward  direction. 
Line  traces  show  the  direction  of  force  which  will  guide  the  local  flow.  Figure  2e  shows  the 
eradication  process  of  the  separation  bubble  formed  due  to  the  flow  pas  t  a  Hat  plate  maintained 
at  +12  deg  AOA.  More  investigation  is  needed  to  ascertain  optimum  electrochemical 
parameters. 


'  0  1  7  1 

arm.*(oti  (e)  Time  evolution  of  separation  control,  near  wall  view 


(d)  Time  averaged  volume  specific  force  field 

Figure  2:  (a-c)  Force  and  velocity  details,  (d)  Time  average  of  streamwise  force  about  the  surface 
of  the  actuator  shows  the  dominance  of  the  forward  (positive)  component,  (e)  Separation  control. 

3.  Polyphase,  multiple  electrodes  and  oilier  control  parameters: 

For  1  mm  thick  Kapton  polyimidc  dielectric  and  12  mm  long  electrodes,  some  effective  control 
parameters  for  active  separation  mitigation  using  asymmetric  dielectric  barrier  discharges  is  studied  by 
considering  the  weakly  ionized  helium  gas  flow  past  a  flat  plate  at  an  angle  of  attack  of  12  deg.  Figures 
3a- f  plot  the  summary  results  after  50  rf  cycles.  It  is  observed  that  the  time  averaged  EHD  force  acts  on 
the  gas  predominantly  downstream,  with  a  transverse  component  towards  the  wall*  The  magnitude  of  the 
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force  depends  on  the  shape  of  the  electrode.  Fig.  3a,  which  also  affects  the  How  actuation.  Fig-  3b. 
Numerical  experiments  [7|  show  that  an  optimum  gap  exists  between  the  actuators  kept  at  a  phase 
difference  for  maxim mn  force  inducement.  Fig.  3c,  and  for  most  effective  separation  mitigation.  Fig.  3d. 
t  he  effectiveness  of  flow  control  becomes  limited  for  more  than  three  actuators  employed  in  phase.  Fig, 


3e-f.  The  effectiveness  can  be  further  increased  if  the  actuators  are  at  an  optimum  phase  difference. 


(a)  EHD  force  distribution  (c)  Force  components  near  electrodes  at  phase  lag 


Stnimwu  dfbiK*.  ijcin| 


(b)  Separation  control  (d)  Actuators  (in  red)  at  optimum  phase  lag  and  distance 


(f)  Flow  along  vertical  lines  4  mm 
downstream  of  the  last  actuator 


(e)  Effect  of  electrodes  on  How  actuation  1  mm  above  wall 


Figure  3:  Effects  of  various  design  parameters. 


4,  Air  Chemistry:  Earlier  we  have  employed  helium  as  working  gas  for  its  simplicity  and  stability  of  the 
discharge.  In  Fig.  4a-d  surface  discharge  has  been  studied  using  air  chemistry  of  formation  of  different 
ion  and  neutral  species  of  nitrogen  and  oxygen  about  an  asymmetric  dielectric  barrier  plasma  actuator. 
Species  with  very  high  recombination  rates  have  been  neglected  for  simplicity.  Response  of  electrons  to 
the  electric  Held  is  very  fast  as  compared  to  ions  due  to  small  mass  and  high  mobility  of  electrons. 
Density  profiles  of  different  species  are  nearly  opposite  to  each  other  at  positive  and  negative  peaks  of  the 
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cycle.  Density  profiles  of  nitrogen  atom  N  and  oxygen  atom  O  are  similar  because  these  do  not  move  due 
to  electric  field.  Density  profiles  of  Nl  and  Oj  are  similar  to  each  other  at  any  specific  time  instance 


because  they  respond  to  electric  field  in  a  similar  manner.  There  is,  however,  difference  in  level  of 
ionization  of  different  species  due  to  difference  in  rate  coefficients,  f  ig.  4c  shows  an  accumulation  of 
charge  very  close  (20  times  Debye  length)  to  the  wall  which  results  in  a  forward  moving  force  vector 
distribution  in  Fig.  4d.  A  detail  air  chemistry  model  with  Poisson-Navjer- Stokes  equations  is  also  under 
development  for  exploration  of  the  effect  of  a  plasma  actuator  operating  in  air  using  realistic  geometry 


and  electrode  configurations. 


(a)  Electron  density  reaches  quasi-steady  state. 


(b)  Distribution  of  negative  oxygen  ion  O'  at  nil. 
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(c)  Evolution  of  q  =  e  (  -  nc  - o  )  at  2GkD 


(d)  Force  vectors  at  the  positive  peak. 


Figure  4:  Solution  details  for  surface  barrier  discharge  using  non-equilibrium  air  chemistry'. 


5.  Integration  of  FDL3D1  and  MIG  for  Flow  Control  Simulation:  Unsteady  END  force  predictions  from 
MIG  have  been  successfully  integrated  into  the  FDL3DI  fluid  code  in  a  collaborative  effort  with 
VAAC/AFRL  researchers.  Data  integration  system  has  been  streamlined  for  smooth  future  transitions. 
Effective  control  of  wing  stall  has  been  investigated  for  a  2  kV  peak-to-peak  voltage  applied  at  5kHz  rf  to 
an  array  of  plasma  actuators  near  the  leading  edge.  Parameters  varied  include  Reynolds  number,  angle  of 
attack,  duty  cycle  frequency  and  actuator  strength  and  location.  For  details  see  Gaitonde  et  ah  |9],  and 
Visbal  et  ah  [II].  Sample  results  plotted  in  Figs.  5a-c  show  that  the  three  dimensional  flow  structure 
obtained  with  a  force  vary  ing  at  radiofrequencies  (5kHz)  is  significantly  different  in  detail  from  that 
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obtained  with  a  steady  actuator.  By  exciting  instabilities  in  the  separated  shear  layer.  One  scale  features 
arise  through  transition  to  turbulence,  which  are  absent  in  steady  actuation. 


Figure  5:  Unsteady  control  of  wing  stall  at  15  deg  AOA  [9,11].  (a)  Plasma  actuators  off.  (b)  Plasma 
actuators  on.  (c)  Mid-span  pressure  coefficients  at  various  plasma  strengths. 

Future  Work 

The  developed  model  is  now  being  employed  for  evaluating  the  flow  control  effect  over  a  vehicle 
forebodv  using  a  polyphase  power  supply  for  various  waveforms  and  duty  cycles.  The  code  is  being 
extended  to  more  complex  air  chemistries  relevant  for  non-equilibrium,  non-therm  a  l  plasma  with 
asymmetric  dielectric  coated  electrode  configuration. 
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Transitions 


MIG  code  has  been  successfully  coupled  with  FDL3D1  code  of  VAAC/AFRL.  MIG  results  for  unsteady 
EHD  force  have  been  incorporated  in  to  the  VAAC  finite  difference  numerical  code  for  simulation  of 
flow  control  over  NACA  airfoil.  Results  were  reported  in  Gaitonde  et  al.  A1AA-2005-5302,  AIAA- 
2006-1205  and  FEDSM2006-98553,  and  in  Visbal  et  al,  AIAA-2006-3230.  Efforts  are  underway  for 
first  principle  simulation  of  vortical  flow  control  about  low  sweep  delta  wing.  AFRL  Contacts  for 
these  transitions  are  Drs.  Datta  Gaitonde  (937-904-4031)  and  Miguel  Visbal  (937-255-2551 ). 
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